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ABSTRACT 


The purpose of this study is to investigate the 
influence of non-environmental and small scale environ- 
mental acoustic factors on the variation of transmission 
loss within areas of the ocean whith exhibit homogeniety of 
large scale environmental factors. The method employed was 
to identify the influencing and non-influencing factors, 
pool over the non-influencing factors if they existed, and 
investigate the nature and magnitude of the effect of the 
influencing factors on the variation of transmission loss. 

The results generally indicated that source/receiver 
depth, the mixed layer, geographic location, and transmit- 
ting range and frequency significantly influenced the 


Variation of transmission loss. 
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I. INTRODUCTION 


In recent years, with the advent of more sophisticated 
SONAR systems, there has been an increasing emphasis placed 
on the study of the acoustic properties of the ocean. To 
be able to predict the capabilities of a particular SONAR 
system, whether it be passive or active in nature, requires 
a complete understanding of the effect of the local acous- 
tic environment on the transmission of sound through the 
water. Both theoretical and observational techniques have 
been utilized in an effort to determine the environmental 
and non-environmental factors which effect the Coameniosron 
Of sound and then to characterize this transmission as a 
function of those factors. 

One study of this nature was the AMOS project performed 
under the sponsorship of the U. S. Navy in the early 1950's 
(ref. 5). In this study, the behavior of sound transmission 
was measured along with local environmental characteristics 
over a wide range of environmental and non-environmental 
factors. On the basis of “hese observations, transmission 
loss prediction equations were developed as a function of 
certain factors such as source and receiver depths, range of 
transmission, mixed layer depth, etc. However, these equa- 
tions are quite insensitive to small scale environmental 
characteristics such as inhomogenieties in the thermal 


structure and irregularities at the air-sea and sea-earth 








interfaces. AS a consequence, transmission losses observed 
operationally are often quite variable with respect to the 
predicted losses on the basis of these equations. 

This study is an extension of the work performed in the 
AMOS experiment in that it investigates the nature of the 
small scale environmental factors which influence the 
observed temporal and spacial variability in acoustic pro- 
pogation. The characterization of this variability as a 
function of these factors will be useful in a variety of 
applications. Some examples are: 

1) Establishing some design specifications 
for new SONAR systems. 

2) Establishing SONAR detection probabilities 
for shipboard use. 

3) Simulation of environmental conditions for 


use with detection simulation models. 








II. A MODEL OF THE ACOUSTIC ENVIRONMENT 


Within sufficiently small areas of the ocean, the large 
scale environmental factors (i.e., thermal, density and 
Salinity profiles) may be considered spacially consistent. 
Such areas may also be considered temporally consistent for 
short durations of time. Since general acoustic behavior 
is a function of these large scale environmental conditions, 
such an area may be considered acoustically homogeneous. 
This, however, does not imply that the acoustic behavior is 
independent of small scale temporal and spacial inhomogeni- 
eties within this large scale environmental structure. 
Hence, the environment may be characterized as a generally 
homogeneous medium which dictates the large scale acoustic 
behavior, but one having small temporal and spacial inhomo- 
genieties which cause variations in acoustic behavior. 

On the basis of this characterization, sound transmission 
loss may be modeled as a stationary stochastic process 
whose mean is a function of the large scale environmental 
factors, and whose autocovariance function is a function of 
the small scale environmental factors. Denote this process, 
where X(t) describes the transmission loss, by X(t);t20. 


It is convenient to characterize this process by its mean 
fe = E(X(t)) 


and by its autocovariance function 








C( A t) = E((K(t) = 7c) (xe) oe 
It is then possible to study the nature of the variability 
which occurs in transmission loss due to small scale 
environmental aspects by considering the environmental and 
non-environmental factors which effect C( At). 

Iwo, particular pabameteus wh Lenwenaisde te ial acme ts 
variability are C(O), the variance of the process and hence 
the amount of variability that would be experienced within 
an acoustically homogeneous region for a given source and 
receiver depths, and At”, the minimum value of At for 
which C(At) = 0, which represents the shortest time inter- 
val between independent observations. These two parameters, 
C(O) and At*”, may be determined for a given process if 
that process is observed either continually or at discrete 
increments of time. Then by observing such processes under 
varying acoustic conditions (such as thermal structure, 
Seographic location, time of day and year, ete.) 1t should 
be possible to determine how C(O) and At* vary as a 
function of the acoustic environment. 

However, it is not necessary to observe the entire 
process to determine C(O) and At™. Speci rien. Ly oie 
the process of transmission loss is observed at two points 
in time separated by Ot and these pairs of points are 
observed many times at varying values of At for the same 


process, or equivalently from processes with the same 
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autocovariance functions, then this information may be used 
to estimate C( At) and hence, C(O) and At*. This may 


be accomplished in the following way: 


Let 


Y( At) = X(t) - X(t+ At) 


represent the difference between two measurements in trans- 
mission loss taken in the same location but separated in 


time by At. 
Then 


Ei Yae Ee = 0 


and 


Var(vy( A t)) E(¥( A t)%) 


Var (X(t) )o+ivar (x +94 2K 


PET Net) 


2 aC EXO) = CAFE) 


The value of Var(Y( A t)) may then be determined for a 
given At from the observations of Y( At) corresponding 
to ee of observations of transmission loss with the same 
autocovariance functions, C( A t). C(O) may then be 


represented by 


C(O) = Var(y( a t*))/2 


me 








where Var (Y( Moa is the value of Var (Y( QA t)) when 
C( A t) approaches zero or when Var(Y(A t)) no longer 
increases with At. At* is then the time interval where 
Enis first occurs. 

If pairs of observations of transmission loss cannot be 
obtained from the same process, then it is necessary to pocl 
information from different processes with the same auto- 
covariance function. The problem here arises in specifying 
the range of conditions under which this may be done. The 
approach taken in this study was to characterize the local 
acoustic environment by its thermal structure where it was 
hypothesized that environments with similar thermal struc- 
tures will produce processes with similar if not identical 
autocovariance functions. This assumes that other variables 
such as transmitter and receiver depths and transmission 


ranges will remain constant. 
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IIl. NATURE OF THE DATA 


The data utilized in this report was made available by 
the U. S. Underwater Sound Laboratory, New London, 
Connecticut, and reflects the raw data collected during 
the Acoustic, Meteorological and Oceanographic Survey (AMOS) 
conducted from June, 1949, through April, 1953. There were 
nine cruises staged during these four years which covered 
the North Atlantic, the Norwegian Sea, and to a lesser 
extent, the Mediterranean Sea. 

Two ships were employed during each cruise where one 
weed as the transmitting platform and =the, other as une 
receiving platform. Each of the cruises was divided into 
several widely spaced stations which served as the focal 
points for data measurement and collection. Within each 
station, the acoustic data was collected at several trans- 
mission ranges between 800 and 30,000 yards. At each 
transmission range the transmission loss measurements were 
observed for four transmitting frequencies at various 
source/receiver depth combinations. The transducer depths 
varied from 20 to 500 feet. The specific data recorded was 
transmission loss, in decibels (db), as a function of 
source/receiver depths and range between ships. Table X 
lists the information that was recorded. 

Accuracy of the data becomes particularly important in 


consideration of the accuracy of any results obtained 


aS 








through the use of that datas  fnestellowing Sic teen com an 
approximate measure of the accuracy of the AMOS data: 
1) time - to the nearest minute 
2) yrange - to the nearest 5 yards 
3) transducer depths - to the nearest foot 
4) water depth - to the nearest 10 fathoms 
5) latitude and longitude - to the nearest minute 
although it is possible that the two ships 
may have drifted as much as a mile or two 
from the recorded position 
6) BT pattern code - reflects the gross thermal 
GOnd1 trons closest EeGgene CimcmeRacmeae 


data was recorded. 
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IV. ASSUMPTIONS 


As indicated in the introduction, this paper deals with 
determining the variation of transmission loss by examining 
the variation of the change in transmission loss as a 
function of the elapsed time between a pair of observations. 
Ideally it would be best if this change in decibel loss 
could be measured from a pair of cbservations with identi- 
cal source and identical receiver depths, and with the same 
horizontal transmitting range, but separated by an interval 
of time, At. Although data of this nature was not available 
from the AMOS experiment, pairs of observations of db loss 
were obtained having the following characteristics: 

1) One of the observed pair measured the 
transmission loss of a signal generated 
at depth d, and received at depth do 
while the other measured transmission loss 
Of a 7signal genera tedwa be dept ime jai 
Be@cmvyed at tAepenr cai. 

2) Identical transmitting ranges between the 
source and receiving ships was not always 


available within a pair of observations. 


Figures (1) and (2) depict the nature of an ideal pair 


of observations and that of the available data. 


The fact that the change in transmission loss was 


measured using transmission paths that ‘crossed’ instead of 
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Available Transmission Paths 


being ‘parallel’ was considered not to be a serious depart- 


ure from the ideal conditions (refs. 


3 and 5). This follows 


from the assumption of a generally constant horizontal 


thermal environment, 


and therefore an acoustically 
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homogeneous environment over the range of operation of the 
two ships on a given station, and theoretical considera- 
tions which postulate that transmission loss over the same 
path from different directions will be the same. It is 
foreseen, however, that this approximation to the ideal case 
of parallel transmission paths will increase the observed 
Var(Yy(At)) for small At due to a contribution caused 
by small scale spacial variability, and may also cause a 
slight underestimation of At”. However, the estimation 
of C(O), the measure of total variation experienced within 
a generally homogeneous acoustic environment, will not be 
effected. 

The second data characteristic of range differences 
between the two observations in a pair was treated in the 
following manner. The observed change in transmission loss 
was adjusted for the effect of this range differential on 
the basis of a piecewise linear approximation of transmis- 
Sion loss vs. horizontal transmission range. This linear 
approximation seemed reasonable in the interest of keeping 
the range adjustment computation simple and yet maintaining 
sufficient accuracy. Plotting transmission loss vs. trans- 
mission range for each source/receiver depth combination at 
each station revealed that the resulting relationships 
exhibited the same general shape and that these relation- 
ships could be approximated by two linear segments with the 


break point at 5,000 yards. 


ING! 








As the slopes of these two linear segments are a function 
of the local environmental conditions, and thus varied from 
station to station, it was necessary to estimate the slope 
for each segment at each station for each source/receiver 
depth combination. This was done on the basis of the trans- 
mission loss observed at the various ranges occupied ona 
given station. Then the change in transmission loss 
observed at a particular station and source/receiver depth 
was adjusted using the corresponding slope of transmission 


loss vs. transmission range in the following manner: 


A db A db -b xX A range 


adjusted zs observed 


where Adb represents the change in transmission loss, )b 
is the slope of the corresponding linear segment, and 

A range is the range differential within a pair of obser- 
vations. 

In some cases insufficient data was available to 
estimate the slope of one or both segments at a given 
station and source/receiver depth combination. In such 
cases the corresponding adjustments were made on the basis 
ofa slope Which was typical of other stations on the same 
and other cruises. Table I depicts a representative cross- 
section of the average slopes from the available data. 
Listed below is a summary of the typical slope values selec- 


ted for use within this study where sufficient data was not 


available for slope calculations. 
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Applicable 


Range Interval Slope Value Cruises 
Above ~O015 Se ec See LON: 
5,000 yards Ligue. 
Below ~0045 7,28; 2, Oe 
5,000 yards ~0100 fel, 
-0280 5 


Admittedly this adjustment procedure is subject to 
error. In particular, the error incurred by using a linear 
approximation in the first place, and secondly in the esti- 
mation of the slope of the transmission loss vs. transmission 
range relationship. It was felt, nowever, that in general 
Pimeetal Error inherent in this procedure could» be 
restricted to at most ‘1.5db by making no adjustments with 
range differentials within pairs of observations of more 
than 500 yards. Errors of such magnitude were considered 
acceptable in light of the original accuracy in measuring 
transmission loss in the AMOS experiment. In actual 
practice, this error was further limited by the fact that 
few adjustments in change in transmission loss were made for 


range differentials greater than 250 yards. 
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V. RESULTS 


A. METHODOLOGY OF INVESTIGATION 

The purpose of this research is to investigate the 
nature and magnitude of the variability in sound propogation 
loss as characterized by C(A t) by investigating the 
Variability in the change in transmission loss as a function 
of time between observations. The calculation of trans- 
mission loss as a function of time, denoted by Y( A t), 
was Obtained from the AMOS data in all instances where a 
pair of transmission loss observations were taken under the 
following conditions: 

1) The compared observations were taken at 
the same station. 

2) The range differential between observations 
was no greater than 500 yards. 

3) The source depth of one transmission equalled 
the receiver depth of the other, and 
vice-versa, thus Hae ree in a compari- 
son of observations with similar but not 


identical transmission paths. 


These calculations of change in transmission loss over 
common paths were indexed by the following environmental 
and non-environmental factors: 

1) Source and receiver depths 


2) Transmission range 
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3) Transmitting frequency 

4) Water depth 

5) Mixed layer depth 

6) Geographical location 

7) Time interval between corresponding observations 


8) Time of day and year 


The investigation of variability in sound transmission 
was confined to the influence of the above factors. To 
determine this influence using the available data, it was 
necessary to pool data collected under varying environ- 
mental conditions. Such pooling is only applicable, however, 
over a range of environmental conditions where the acoustic 
variability is characterized by a common autocovariance 
function. Hence, it is necessary to determine how to iden- 
tify the range of conditions which lead to a common auto- 
Sevarlance function. In particular, it is necessary to 
determine which factors, both environmental and non- 
environmental, influence the nature and magnitude of the 
variation of propogation loss in the ocean and then pool 
Over the entire range of those factors which did not so 
influence this variation. 

The procedure followed was to plot Var(y( At)) vs. 

At using the entire sample set of data and thus represent- 
ing the entire range of environmental and non-environmental 


conditions encountered. This entire sample of data was 
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then subdivided into sub-samples on the basis of the 
source/receiver depth combination and Var(Y(.A t)) vs. 
At was plotted for each sub-sample. By this procedure, 
it was determined that the nature and magnitude of 
Var(Y( A t)) depended upon the sub-sample (i.e., source/ 
receiver depth combination) considered. Hence, it was 
deemed inappropriate to pool the data over source/receiver 
depth combinations when considering the influence of other 
factors on Var(Y( A t)). 

To determine if the thermal structure influenced 
Var(Y( A t)), each of the above sub-samples was further 
subdivided on the basis of the location of the receiver and 
source depths with respect to mixed layer depth. On the 
basis of this subdivision, it was determined that the posi- 
tioning of the source and the receiver with respect to the 
mixed layer depth also influenced the nature and magnitude 
mee Var(Y( A t)). 

The influence of geography and transmission range were 
also investigated in the same manner by subdividing the 
sample data on the basis of these factors. As it was 
necessary to keep the respective sub-samples as large as 
Bessible when considering a given factor, it was found that 
certain sub-samples should not be further subdivided. Such 
was the case for some of the less representative source/ 


receiver depths. 


Ze 








By proceeding in this manner, it was determined that, of 
the factors being considered, those which were judged to 
have a significant influence on the nature and magnitude of 
Var(y( A t)) were: 1) source/receiver depth, 2) location 
of the source and receiver with respect to the mixed layer, 
3) geographic location, and 4) transmission range. It was 
also determined that the transmitting frequency had a marked 
effect on the magnitude of Var(Y( A t)). The actual nature 
of these influences will be discussed later. 

It should be emphasized that Var(Y( At)) is twice 
the value of C(O) and thus Var(Y(t)), the variation os 


transmission loss. 


B. LIMITATIONS EFFECTING RESULTS 

There are two considerations which tended to limit the 
results of this paper. The first was that the amount of 
data available to estimate Var(Yy( At)) as a function of 
At was limited in some of the finer subsamples defined in 
the previous section. Hence it was round necessary, in 
order to obtain a reasonably precise estimate of 
Var(Y¥( A t)), to pool the available data into several con- 
secutive 15 minute intervals. 

The second consideration was that, due to physical 
limitations imposed on observation of the data during the 
AMOS experiment for a given source/receiver depth combina- 


tion, the majority of the available data was restricted to 
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at most three 15 minute intervals. These three intervals, 
however, were not necessarily the same for each source/ 
receiver depth combination, but tended to the small values 
of At for closely spaced source and receiver depths and 
to the large values for widely spaced source and receiver 
depths. 

The implication of these data limitations was that A t* 
could not be determined with a precision greater than 15 
Minutes. However, this was thought not to be a serious 
limitation. A possibly more serious limitation is that, 
because of the dispersion of the data as a function of At, 
the value of Var(y( At)) cannot be estimated with equal 
precision from interval-to-interval and from depth 
combination-to-depth combination. In particular, fora 
given source/receiver depth combination, there may be little 
pamno data available £or the interval when Var{(y(A £)) 
reaches its maximum value. In general it is felt that this 
did not happen, although there is no way of verifying this 


at the present time. 


C. FINDINGS 
1. Distribution of Transmission Loss 
As suggested in the introduction, knowledge of the 
nature of the variability in transmission loss may be 


utilized to produce a realistic simulation of the acoustic 
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environment. To do this it is necessary to know not only 
the magnitude of this variation, but also its distribution. 

The typical assumption that this distribution is normal 
is evaluated on the basis of the available data for a 
representative cross-section of subsamples of the data. 
The basis for this examination is the Kolomogorov-Smirnov 
one-sample goodness-of-fit test (ref. 2). Table II gives 
the results for all subsamples considered, and the range of 
conditions which define the tested subsample. Also given 
ae D4, the maximum absolute aqifference between sample 
cumulative distribution and a normal cumulative distribution 
With mean zero and sample variance of the subsample. The 
critical value has the same connotation usually associated 
with hypothesis testing and represents the level of signifi- 
cance at which the null hypothesis of no departure from 
normality would be rejected. 

The general conclusion on the basis of these results is 
taet ue aera distribution is an adequate approximation 
for the distribution of transmission loss. 


2. Factors Effecting the Variation of Transmission Loss 


a. Source/Receiver Depth 
The first attempt at determining which environ- 
mental and non-environmental factors influenced the nature 
and magnitude of Var(Y( At)) was in the area of 
source/receiver depths. The available data was subdivided 


on the basis of several source/receiver depth combinations. 
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The results thus obtained proved to be quite noisy and it 
was decided to employ the same procedure for cruise 10 only 
Since the thermal structure appeared more consistent on that 
cruise. The results of this procedure seemed satisfactory 
for exposition purposes and are summarized in Table III. 

It is interesting to note that not only does Var(Yy(At)) 
increase as time increases, but that this variation also 
increases with an increasing depth differential between 
source and receiver. Also noted was that Var(yYy(A t)) 
tended to decrease as the source depth increased. The 
latter is probably attributable to surface effects. 

b. Mixed Layer Depth 

Once it had been decided that the different 
source/receiver depth combinations each influenced 
Var(Yy( At)) differently, the next step was to consider 
the position of the source and receiver with respect to the 
depth of the mixed layer. This was accomplished by further 
subdividing the data into two subsamples on the basis of 
whether both the source and receiver were in the mixed 
layer, or at least one was below the mixed layer. The 
rational behind this subdivision was to determine if the 
variability of transmission loss was different when the 
transmission path was above the thermocline or across the 
Enermoc line. 

Table IV shows the results of this investiga- 


tion and clearly indicates that Var(y(At)) is 
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considerably less within the mixed layer. The differences 
observed may in part be a function the existence of internal 
waves. 
Cc. Geographic Location 
Geographic location was found to have a 

Significant effect on the variation of transmission loss, 
To facilitate investigating the nature of this effect, areas 
were selected which were thought to have approximately the 
same large scale environmental characteristics. The areas 
chosen were: 

1) 50-80 N, all lor.gitudes 

2) O-50 N, O-60 W 


3) O-50°N, 60-90 Wi 


Four source/receiver depth combinations were then 
picked which would ensure that sufficient pairs of observa- 
tions would be available within each of these depth combina- 
tions. The subsamples have been further subdivided accord~ 
ing = the position of the source and receiver with respect 
to the mixed layer. Table V depicts Var(Y¥( A t)) for the 
geographic areas and the four depth combinations selected. 
The variation of transmission loss for each combination 
pooled over all latitudes and longitudes has also been 
included in the table for comparative purposes. 

As indicated by the entries in Table V, there 
does not appear to be any consistently predictable direc- 
tion or magnitude associated with Var(Y( At)) across the 
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three geographic locations. However, the variations do seem 
meechange significantly both from the "all lLatitude- 
longitude" case, and across geographic areas. This would 
seem to substantiate the hypothesis that geographic loca- 
tion is a significant factor effecting Var(y( At)). The 
data in Table VI, which breaks the variation of trans- 
mission loss down into 15 minute intervals, also appears 
to substantiate this contention. 

d. Transmitting Frequency 

The effect of transmitting frequency on 

Var(y( 4 t)) ais apparent; increasing the frequency in- 
creased the magnitude of the variation. Two frequencies 
were selected to investigate this hypothesis, 8 and 25 KC. 
Table VII is a summary of Var(y( At)) for both frequen- 
cies for several source/receiver depth-geographic location- 
mixed layer depth combinations. The increase in Var(Y(At)) 
is evident. 

e. Transmission Range 

The transmission range appeared to have a 

Significant effect on Var(Y¥( At)). Table VIII shows this 
variation ror range intervals of O-10 Kyds and 10-30 kKyds 
for pairs of observations at the 8 KC transmitting frequency 
pooled appropriately according to source/receiver depth, 
geographic location and mixed layer depth. An examination 
of this data suggests that Var(y( A t)) was generally 
lower for greater ranges and that estimates of variation 
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in the mixed layer was less consistent than for those 


transmissions crossing the layer. 


Se Investigation of Large Values of Change in 
Transmission Loss 


Because of the limited amount of data in some 
subsamples considered, the estimator of Var(Y( At)) was 
quite sensitive to a not uncommon occurrence of dispropor- 
tionately large values of the change in transmission loss. 
In an attempt to determine the nature of the circumstances 
Which lead to such observations, a sample of observed 
Changes in transmission loss whose absolute value exceeded 
9db were examined as to the conditions under which they were 
observed. This sample consisted of all instances for 
20/100, 50/100, 50/250 and 100/500 foot source/receiver 
depths at 8 KC transmitting frequency, and 20/100 and 100/500 
foot source/receiver depths at 25 KC transmitting frequency. 
Table IX is a summary of these data. 

On examination of these data it is concluded that, 
except for a marked effect of frequency and a suggested 
effect of time of day, no common single or set of circum- 
stances is related to the large changes in transmission 
loss. This does not mean that such a set of circumstances 
does not exist, but that it is not apparent what such a set 
would be based on in consideration of these data and this 


specific set of indices. 
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This latter point has lead to the consideration 
that there may be a more appropriate set of indices with 
which a measured change in transmission loss may be classi- 
fied. The nature of such a set of indices is suggested by 
considering the results of the present study. That is, that 
Var(Y¥( A t)) varied with the proximity of the source or the 
receiver to the surface, the depth differential between the 
source and receiver, the relationship O£ the source and 
meeiver to the mixed layer, geographic Location, etc. 

These factors are exactly those which influence the mode ox 
Peansmission of a signal in the general transmission loss 
equations obtained in the AMOS study (ref. 5). 

In that study it was determined that transmission 
loss depended on the mode of transmission. In particular 
On whether transmission was by 

1) Direct path 

2) Single surface bounce 

3) Multiple surface bounce 

4) Leakage by diffraction or scattering 
5) Depressed sound channel 

6) Bottom bounce 

On the basis of the results of the present study, 
it appears that the nature and magnitude of the variability 
in transmission loss might also vary with transmission mode. 


Further, it is envisioned that subdividing the available 
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data on the basis of transmission mode will yield a more 
uniform distribution of At for a given subsample and, in 
general, increase the amount of data available for the esti- 


Mation of Var(Y¥( A t)) for any particular interval of At. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The aim of this study was to determine the environmental 
and non-environmental factors which influence the nature 
and magnitude of the variability in transmission loss, and 
then to characterize this variability as a function of these 
Influencing factors. 

On the basis of this study it was found that the 
variation in transmission loss was influenced by the depth 
of the source and receiver of the transmitted signal, the 
relationship of the source and receiver depths with the 
mixed layer depth, geographical location, and transmission 
range. 

The nature of this influence was that the magnitude of 
variability increased with increasing transmitting frequency, 
with increasing depth differential between the source and 
receiver, with the increasing closeness of the source and 
receiver to the air-sea interface, with decreasing trans- 
mission range, and when transmitting across the thermocline 
as opposed to transmitting exclusively within the mixed 
layer. 

Te time interval between uncorrelated observations was 
found to be shortest near the surface as opposed to at 
depth, to increase with increasing source/receiver depth 
differential, and to decrease when transmitting across the 


thermocline as opposed to transmitting entirely within the 


mixed layer. 
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The magnitude of the variability in transmission loss 


experienced under the varying conditions ranges over one 


2 in the near horizontal trans- 


2 


order of magnitude from 10 db 
mission in the mixed layer at 8 KC to 100 db” crossing the 
Emermochine in the 25 KC case. Such figures sare onlyeecer 
mates and as such are subject to statistical variation, but 
the trends in magnitude of the variability appear to be 
consistent in the various cases considered. Therefore, the 
magnitudes of the variability determined in this study may 
be used as approximations of the true variability which will 
be experienced in the corresponding environmental situation. 

In conclusion, it is observed that the factors which 
influenced the nature and magnitude of the variability in 
transmission loss are exactly the same as those found in the 
AMOS study to influence the mean transmission loss (ref. 5}. 
imeparticular, the factors of source and receiver dépen, 
mixed layer depth, transmission range and surface tempera- 
ture were used to determine the mode of transmission of a 
given transmitted signal. The transmission loss expected 
was then calculated using an equation peculiar to this speci- 
fic transmission mode. 

In a like manner it was suggested that the nature of the 
variability in transmission loss may be more adequately 
classified and the magnitude of the variability more pre- 


cisely specified if the observed change in transmission 


losses were further indexed and examined by transmission 


a5 








mode. It is therefore recommended that this be the next 
area Of investigation in the study of the nature and magni- 


tude of the variability of transmission loss in the ocean. 
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APPENDIX A 


I. DESCRIPTION OF DATA 


A. TAPE RECONSTRUCTION 
The information as received from the U. S. Underwater 

Sound Laboratory was in the form of 9, seven-track tapes 
Containing appropriate spacial, temporal, and acoustic 
information. These tapes have now been consolidated into 
two nine-track tapes suitable for mounting on the standard 
IBM tape drive units. One tape contains only the bathy- 
thermograph (BT) information while the other contains all 
Other necessary data including acoustic measurements 
Meg@teanent to this report. Table X displays the information 
that is recorded on this latter tape. The following informa- 
tion has been omitted from each record on the revised 
meoustic tape, but is still available on the original seven-— 
track tapes: 

1) Ambient noise indicator 

2) Wind force (beaufort) 

3) Unknown contents following water depth entry 

The acoustic tape contains data for cruises 5, 7, 8&, 9, 

io ll, and 12. Each record is 7/2 characters in length 
although each only utilizes 60 character spaces. The addi- 
tional 12 spaces have been left blank to allow room for any 


additional information deemed appropriate and to facilitate 
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placing this data on a user disc if desired. No special 
records have been added to indicate the end of a cruise or 
station. It should be noted that all field entries are in 
terms of integer numbers with the exception of fields 5, 
24, and 27 which will require character format (A~format 
under FORTRAN IV). 

The data for cruises 2 and 3 were deemed inappropriate 
and therefore not included on the acoustic tape. The 
reason was that the source/receiver depths and two of the 
four transmitting frequencies do not correspond to the 
requirements of this study. That is, both cruises show 
source depths set at a constant 15 feet while the receiver 
depth is varied. Consequently, none of this data would be 
accepted by the computer program. The program requires the 
Comparison of one record with source depth at d, feet and 
mie receiver depth of dj feet, with another record having 


the source at depth dy and receiver at depth dj. 


B. COMBINING ACOUSTIC DATA 

An examination of the original tapes revealed that 
either 1, 2, or 3 acoustic measurements were taken for each 
Maztacular location, time, and range. Hence, some method 
must be employed to combine these sets of data in order that 
they be consistent. One approach would have been to take an 
average each time a grouping of two or three measurements 


was encountered. Another approach, and the one subsequently 
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adopted, would be to simply sum the acoustic data for 
groupings of three measurements, sum and multiply by 3/2 
for groupings of two, and triple any single entries. fMThe 
resulting figures would then be entered on the revised 
acoustic tape in place of the old transmitting frequency 
data, giving a single entry for each particular location, 
mame, and range. This latter,.method was utilized in the 
interest of conserving valuable computer time, and in 
Consideration of computational and programming ease. In 
Summary, all acoustic field entries at each of the four 
Pmeansimitting Erequencies must be Givided by 3 before ehey 


Can be used in any computations. 


C. BATHYTHERMOGRAPH PATTERN CODE AND ITS USE 

The bathythermograph (BT) pattern code is a two digit 
number and will be found in fields 18 and 20 on the acoustic 
ape . 

Interpretation of the code has been simplified for use 
weemin this study. Figure (3) indicates the general nature 
of the thermal structure. To use the figure, enter with 
mae first digit of the code across the top. This gives a 
picture of the approximate temperature gradient down to 100 
feet (associated depths may be read in the ‘lst digit’ 
column). Next, enter with the second digit. This gives a 


picture of the approximate temperature gradient from 100 to 
250 feet. Connecting the two pictures thus obtained results 
in a complete recap of the thermal conditions from 0 to 250 
feet. 

37 





Depths (ft) Pattern Code Digit 


et react Pyoval (Glatlefshe i 2 3 os 2 6 de 
0 100 
30 150 
50 200 
100 250 


Figure 3. —- BTePattern Code Interpretation 


This BT pattern code information was used to establish 


the position of the mixed layer. 
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APPENDIX B 


I. EXPLANATION OF COMPUTER PROGRAM 


This section has been included in the report in order 
that the program may be used as a base for further studies 
in the field of temporal and spacial effects on the varia- 
tion of decibel loss underwater. The entire program has 


been written and compiled under FORTRAN IV, G-level. 


A. GENERAL EXPLANATION 

The program is designed to read one data record at a 
time from the acoustic tape, check the appropriate data 
against the parameters established by the input array, and, 
if possible, pass the entire record to subroutine COMP1. 
Subroutine COMP1 then separates the records passed to it 
into four categories according to range between the trans- 
ieecang and receiving ships, the range break poine, and 
source and receiver depths. For example, if the source/ 
receiver depths are either 20/100 or 100/20 feet and the 


range break point is 5,000 yards, the categories are: 


Source/Receiver 
mca cegory, Depth Range Interval 
I 20/100 = 5,000 yards 
II 20/100 >5,000 yards 
en 100/20 = 5,000 yards 
ay 100/20 > 5,000 yards 
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Once this data has been collected for a single station, the 
appropriate slope information is computed by subroutine 
RGRESS and the next station is then processed through the 
Main program, COMP1, and RGRESS. 

At the end of each cruise, the processed records are 
sent to subroutine COMP. This routine computes the change 
in transmission (db) loss in a sequential manner. That is, 
a record can never be used in a computation more than once. 
The computation of Adb loss proceeds in the following 
manner: 

1) At any given time, two adjacent records are 
“considered. 

2) The program computes the difference between 
ranges of the two records. If this difference 
is greater than the specified amount 'Z‘', then 
no computation is made and the next data set 
is picked up. 

3) If the range difference is less than or equal 
to the specified amount, Z, then the change 


in db loss is calculated. 


of the. two records being considered at any particular time, 
the decibel loss of the second is adjusted in the direction 
@f the first record with respect to the difference in range 
associated with the two records. This is the point at which 


the slope information passed from subroutine COMP1 is util- 


ized. The only exception to this computational procedure 
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is when one range falls below the break point and the other 


above. In that event, both db losses are adjusted to the 
range break point. Once a computation has been completed, 
the change in decibel loss and the time difference between 
the two records are stored in appropriate arrays. At the 
end of each cruise, this information may be plotted at the 
users option. Also within this subroutine are two arrays 
which collect the same information for all cruises to 
facilitate a cumulative plot. 

At the end of the tape, the program calls subroutines 
SIGMA and BTLAYR. SIGMA simply sorts the Adb losses with 
respect to their associated time increments into 5 and then 
15 minute intervals, computes the variance within each 
interval, and plots the variance vs. the change in time. 

It also calculates the variance and standard deviation over 
all data passed to it at the end of each tape iteration. 


BTLAYR is designed to consider all comparisons which 
lead to computation of Adb loss in subroutine COMP, and 
separate the pairs of observations into one of three cate- 


gories: 1) above the layer, 2) transmitter and receiver 
both below the layer, and 3) all transmissions with the 
receiver and source neither above or below the thermocline. 
Once all records have been categorized, the change in db 


loss vs. change in time may be plotted at the users option. 


Next subroutine SIGMA is called, which computes the 
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appropriate output as described above for each of the three 
transmission categories. 

Subroutine SIME solves simultaneous linear ee 
necessary to subroutine RGRESS. OSPLOT is a standard IBM 
plotting package and is used exclusively for plotting within 


the program. 


B. INPUT SPECIFICATIONS AND DESCRIPTION 

An input array has been established in order to lend a 
Sertain amount Of flexibility in filtering data through 
specific parameters. This information must be punched on 
computer cards for inclusion at the end of the program and 
Will be either integer or alphabetic in nature. There must 
be six entries on each card with the last letter of the 
numeral entered successively in columns 6, 12, 18, 24, 30, 
36, 42, 48, 54, 60, and 66. The last 9 entries on the last 
Card may be left blank. Table XI fully explains what 
information is required and gives a general description of 
the input variables. The input data is read into the NK 
array. Note that NK(15) and NK(17) require a character 
format. Tables XII and XIII describe the arrays and vari- 


ables used within the computer program. 


C. READING DATA FROM THE ACOUSTIC TAPE 
The acoustic tape data format is described by Table X. 


There is one peculiarity which exists within the formatting 


structure. Note that the range between the two ships is 
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divided into two separate field entries. In this case it 
is necessary to compute the range in the following manner 


(tape records are loaded into the IN array): 
Range = 1N(9) + 1,000 X IN(8) 


All data on the tape is integer in nature except fields 5, 


6, 24 and 27, which are alphameric. 


D. ESTABLISHING THE MIXED LAYER DEPTH 

Within each pair of compared data sets are four BT 
pattern codes, one per ship per record. Since there exists 
a Change in time and a range differential between data sets, 
it seemed appropriate to combine the depths associated with 
the pattern codes in some manner. One acceptable method and 
the one adopted, takes an average of the associated thermal 
layer depths across records first for the transmitting ship, 
and then the receiving ship. This averaging is done in the 
Same manner for determining both the top and the bottom of 
the thermocline. Once this has been accomplished, the mini- 
mum of the two top of the i AER and the maximum of the two 
bottom of the layer averages are selected to represent the 
upper and lower bounds of the thermocline respectively. 

As an example, consider some BT pattern code XjX9, 
where it was desired to determine both the top and the 
bottom of the thermocline. Computing the top and bottom of 


the thermocline was done in the following manner: 
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1) Determining the top of the thermocline (enter 





with Xj) 
Xy Depth (ft) XQ Depth (ft) 
1 O 1 100 
2 O 2 100 
3 see X5 3 100 
ot O & 100 
5 30 5) 150 
6 30 6 150 
i 50 7 200 
8 see Xp 8 600 


2) Determining the bottom of the thermoline 


(enter with X5) 


X95 Depth (ft) XL Depth (ft) 
1 300 1 100 
2 200 2 50 
3 150 3 30 
4 300 4 100 
5 300 » 50 
6 300 6 100 
7 300 7 100 
8 see Xj 8 600 


Note that when checking for the top of the thermocline that 


when X, = 3 or 8, this implies looking at the second digit. 


Aa, 








Likewise, when determining the bottom of the thermocline 


and X5 = 8, then look at the first digit. 
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TABLE I. 


Average Slopes Pooled Over All Possible Data 


Slope 
* 
di /d d/a,* 
Source/Receiver Cruise 

Depth (ft.) No. £5 Kyds 75 Kkyds *5 kyds 75 Kyds 
50/200 5 -0O03I13m)) .06622 )ecoOlcge 
7 .001978 

8 moles meen 1 

9 .004476 .001040 

10 .003580 .001894 

aL .010430 .000683 
30/150 9 -004950 .000759 .000358 
10 .004975 .001400 .004202 .001356 
20/100 5 .028020 .002022 .023260 .003018 
y 20019581 002236 
8 .005097 .001348 .004221 .015560 
9 .009709 .001062 .008700 .001550 
iL .009000 .001000 #.010100 .001000 
50/400 8 .003831 .001400 .003430 .001358 

9 .004455 .001191 

10 .003607 .001924 

11 Aonlibalal(ey. * stelererse7 

100/400 5 .000599 

| 7 .000930 
8 .004718 .001453 .001964 .001313 

9 {006560 me OOM 29 

10 .003980 .001434 

id .016850 .001105 


*source/receiver depths 
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TABLE iV: 


The Effect of Mixed Layer Vone Var "Gen t)) 


Source/Receiver 
Depth (ft.) Above™ 
20/100 14.8 
20/500 D7 se 
30/150 BOS han 
50/100 2055 
507250 Sy 
50/400 LSi.10 
50/500 14.6 
100/400 a ee, 
100/500 12.0 


ee 
Crossing 


28.4 


5Ge/ 


24.9 
256! 
Zo) 
34.3 
25.4 


41.5 


# . : 
Both source and receiver in the mixed layer 


AE ; ; 
At least source or receiver below the mixed layer, i.e., 


crossing the thermocline 
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TABLE V. 


The Gross Effect of Geographic Location 


on Var ( A db loss) - 8 Kc 


Source/Receiver Latitude Longitude > 
Depth (ft.) (degs) (degs) Crossing Above a 1 
20/100 eal abl 26030 (ae) > allteee (ons) 
50-80N All 23242 (amie) 2.60(2) 
0-50N 0-60W 37 1G0(G5)n menos Olle) 
O-50N 60-90W 8.30(5) Se kot 2) 
50/100 All All 24.95(216) 26253 (067) 
50-80N All 20200 (155) 6. 77X55) 
O-50N O-60W 29.96(79) Zia Solas) 
O-50N 60-90W 25.95 (i) 29.10(29) 
50/250 All All 25/1 GaSe) 7.53 (48) 
50-S0N All 23554 (94) 4.86 (39) 
O-50N 0-60W 20m 5213 5) 24.61(6) 
O-50N 60-90W AS 2622) 7.89(3) 
100/500 All All 41.58(116) 11.97 (40) 
50-80N All 39.39 (87) 2.51) (345 
O-50N O-60W 
O-50N 60-90W 49.87 (28) -43 (3) 
*{ ) — Indicates sample size. 


* 


*at least source or receiver below mixed layer. 


***poth source and receiver in the mixed layer. 
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Field 





WOON OA UPWNH-E 


Columns 


AVN SPA D.C. 


Acoustic Data Records 


HHH 


» 


HHHHHHHHHHHHRARARRRRARARRaAa aR 
BEONPNNNEFNEFWWWWWWHEWYHNNDND 


Format 


Hm NO DO 


NO 


Description 


Cruise Number 

Station Number 

Hours (GCT) 

Minutes (GCT) 

Sign Of @imers4Zone 

Time Zone 

Day 

Month 

Year 

Range (thousands of yards) 
Range (hundreds of yards) 
Receiving Ship 

South depth (ft.) 

Receiver depth (ft.) 
Propogation loss (db) 2.2 KC 
Propogation loss (db) 8 Kc 
Propogation loss (db) 16 Kc 
Propogation loss (db) 25 Kc 
Sea State ) 

BT Pattern Code ) Ship 1 
Sea State ) 
BT Pattern Code ) 
Degrees (latitude) 
Minutes (latitude) 
N or S 

Degrees (longitude) 
Minutes (longitude) 
E or W 

Water depth (fathoms) 


Sem yey 


64 


10 
eA 


a2 
1.3 


14 


TABLE XI. 
Input Array Description - NK(I) 


Description 


Cruise number, Zero implies all cruises, otherwise 
specify cruise number desired. 


Station number. Zero specifies all stations, 
otherwise specify station number desired. 


Time (GCT). Zero implies all times, otherwise 
specify lower bound on time interval desired. 


Time (GCT). Same as NK(3) except specify upper bound 
On desired time interval. 


Time zone. Zero implies all time zones, otherwise 
specify desired time zone. If negative time Zone, 
fill in first with negative sign, otherwise leave 
first character blank. 


Month. Zero implies all months, otherwise fill in 
lower bound no month interval desired. 


Month. Same as NK(6) except specify upper bound on 
month interval desired. 


Year. Zero implies all years, otherwise specify 
year desired. 


Range. Zero implies all ranges, otherwise specify 
lower bound on range interval. 


Range, Same as NK(9) except specify upper bound. 


Receiving ship. Zero implies both ships, 1 implies 
Ship No. i, 2 implies Ship Woe. 


Source depth. Specify source depth desired. 
Receiver depth. Specify receiver depth desired. 
Latitude. Zero implies all latitudes. This is a 
packed word which expresses latitude in degrees 


and thousandths of degrees. For example: 40.817 
degrees will be put in the array as 40817, and 
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(eont' d) 


iS 
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al 
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ite, 


20 


ark. 


Ze 


218 


24 


25 


26 


TABLE XI. (Gomme) 


31.7 as 31700. This will always be the ‘lower 
limit. If working with a N/N or S/S situation, 
enter the smaller of the two latitudes. If 
working with a N/S situation, enter the north 
latitude. 


Nor S. This will always be an N for the N/S 
Sltvation. 


Latitude. Zero implies all latitudes. This is 
a packed word (see NK(14)). This will always be 
the upper bound on the latitude interval. If 
working with a N/N or /S/7S situatwonm enter sene 
larger of the two latitude limits. If a N/S 
situation, enter the south latitude. 


Nor S. This will always be S for N/S situation. 


Longitude. Zero implies all longitudes. This 
is a packed word (see NK(14)). This will always 
be the lower limit. If working with E/E or W/W 
Situation, enter the smaller of the two longi- 
tudes. If E/W situation, enter the east 

long. Eudes 


E or W. This will always be E for E/W situation. 
Longitude. Zero implies all longitudes, This is 
a packed word (see NK(14)). This is always the 
upper limit on longitude] | bneer wehew [amc emer 
the two longitudes here if a E/E or W/W situation. 
Enter west longitude for E/W situation. 

Eor W. If E/W situation, enter W here. 


Water depth. -l implies all depths, otherwise 
specify lower limit. 


Water depth. Same as NK(22) except specify upper 
MEIER Te te 6 


Transmitter frequency. 13, 14, 15 and 16 imply 
2.2, &, 16, and 25 ke tespcerum cry. 


Number of stations in cruise 5. 


Range differential between pairs of observations. 
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ZS 
Zz 
30 
Sil: 
3 yd 
33 
34 


35 


TABLE XI. (Cont'd) 


Number of stations cruise 7. 

Number of stations cruise 8. 

Number of stations cruise 9. 

Number of stations cruise 10. 

Number of stations cruise ll. 

Number of stations cruise 12. 

Enter first cruise number to be used. 

Specify. range brea eewmr. 

Number of repetitions of program desired. Must 


have a complete set of input data cards for each 
LECratJoOn. 
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Arrays 


A 


Al 


adh 


AU 


AV 


Bl 


BT 
BU 

BV 

DB1 
DB2 
DB3 
DB4 
DBO 


DB 


TABLE XII. 
Description Of Amma a 


Description 


Stores the time difference twithin, pareomor 
observations by cruises. 


Accumulates time differences over all cruises. 


Stores time differences for pairs of observa- 
tions that have both receiver and transmitter 
neither above or below the thermocline. 


Stores time differences for pairs of observa- 
tions that have both receiver and transmitter 
above the thermocline. 


Stores time differences for pairs of observa- 
tions with both transmitter and receiver below 
the thermocline. 


Stores the difference between adjusted trans- 
mission losses within pairs of observations 
by cruises. 


Accumulates adjusted transmission loss 
differences over all cruises. 


Same as AT for transmission loss differences, 
Same as AU for transmission loss differences. 
Same as AV for transmission loss differences. 


* 
Decibel loss for Case I. 


Decibel loss for Case ac 
Decibel loss for Case III. 
Decibel loss for Case Iv. 
Combines db losses for Cases III and lis” 


Combines db losses for Cases I and ii © 
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IK 


IN 


LUX1L 


LLX2 


LUX3 


LUX4 


NK 


PDB 
PRNG 
RNGL 


RNG2 


TABLE. XiGs (Gont > d) 
stores first and secondsdiot¢s fovea t pace mn 
code. 


Array MK is loaded into this array one record 
at a time for computations in subroutine COMP. 


Loaded sequentially with records from the 
acoustic data tape in the main program. 


Establishes midpoints of five-minute intervals 
in subroutine SIGMA. 


stores BY pattern code for Shipenereeenee are 
data set in pair. 


Stores BY pattern code for "ship ele wean so. 
data set in pair. 


Stores BT pattern code for Ship No. 1 and 2nd 
data set in pair. 


Stores BT pattern code for Ship No. 2 and 2nd 
data set in pair. 


Stores records from acoustic tape that have 
the same cruise no. 


Temporary storage for acoustic tape records 
when a different cruise number other than the 
one being considered is encountered. Will be 
placed in the MK array after computations are 
completed for the previous cruise. 


Counts entries in each five-minute interval in 
subroutine SIGMA. 


The input parameters. 


Counts entries in each 15-minute interval in 
subroutine SIGMA, 


Stores DBO for plotting at end of each cruise. 
Stores RNGO for plotting at end of cach Cruise. 
Stores ranges for Case Tae 


x 
Stores ranges for Case II. 
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RNG3 


RNG4 


RNGO 

SLOPE 
SLOPE 
SLOPE 


SLOPE 


XRNG1L 
XRNG2 


XRNG3 


XRNG4 


XRNGO 


XRNG 


XS 


*Refer diagram on Page 39. 


Stores 


Stores 


Stores 


Stores 


Stores 


Stores 


Stores 


Stores 


Stores 


Stores 


Used 
Used 
Used 
Used 
Used 
Used 


Used 


41591 


abel 


sgl 


alg: 


in 


In 


TABLE Xia 


ranges 
ranges 
ranges 
ranges 
ranges 
slopes 
slopes 
slopes 


Var (Y ( 


(Cont'd) 


for Case III.* 

for Case. mee 

for Cases I and ey di 

for Cases ule. and an” 

for Cases III and Athy 

for Case II” by Station numoer. 

for Case III* by station number. 
for Case Iv* py station number. 


t)) within five-minute intervals. 


top and bottom Of “Ehernec line 
information for pairs of observations. 


PlLOteiIngd RNG 


plotting RNG2. 


PLOLEING =RNGSs- 


plotting RNG4. 


plotting RNGO. 


plotting RNG 


in computing Var(y(At)) in subroutine 
SIGMA. . 


Var (Y(At)) 


in subroutine 


Within each 15-minute interval 
SIGMA. 
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Variable 


TA 


IJK 


IQ 


JJK 


JIL 


be 


LL 


LUL 


TABEE <2 in. 


Description of Variables 


Description 


Set equal to the station number. Used to 
check for the end of a station. 


set to the cruise number under consideration. 
USCC anieasincasnic migige 


Counter for arrays collecting sda rosevomemeaan 
cruise, 


May be 1 or zero. When set to zero by the 
program, 16 Stops ENG 1teraGion dee ene seme son 
the cruise specified by NK(l). 


A switching mechanism that allows computations 
to proceed for cruise 12 .aeeenese ieee = een 
aCOUSeUC Udo ca mite. 


Specifies the number of Iterations within scsen 
run of the program. 


Cruse Numpeias 
Sets the transmitting frequency. 


Specifies the number O£ Stations inva 
particular cruise. 


Allows for zeroing certain variables in sub- 
routine COMP, at the beginning of each 
Be ie terelite niece tals 


The BT pattern code for Ship ff ama 1st set 
in data pair. 


The BT pattern code for Ship 2 and lst set 
in date pare: 


The BT pattern code for Ship 2 and Ist set in 
Sleliesl joveiilie 


The BT pattern code for Ship 1 and 2nd set 
in data padre 


ab 





Variable 


LX4 


MM 


SIG 


N KK * 5 


TABLE cn in (Conta) 


Description 


The BT pattern code for Ship 2 and 2nd set 
in data pair. 


Counter for the MK array. 
Counter for the MMK array. 


Counter for arrays {collecting wed rage. 
a Meee nuisescr, 


Set to the cruise number under consideraion. 
Used to check for the end of a cruise. 


Sets the station number for use with the 
SLOPE arrays. 


Standard deviation of Y(t) 
VaEety( ty 

Source depth 

Receiver depth 


Maximum allowable range differences between 
pairs of observations. 
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